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Currently available investigations of the isotope effect in Fe-based superconductors are all based on 
macroscopic resistivity and susceptibility measurements and results are highly controversial. Here we 
present an investigation of the Fe isotope effect in the (Ba, K)Fe 2 As 2 superconductor by 
temperature-dependent EXAFS in the framework of Einstein model. This method separates the static 
structural isotope contribution and the local lattice isotope contribution. Results point out the isotope 
substitution has a negligible effect on inter-atomic distances. On the contrary, both the local vibrations of the 
Fe-As and the Fe-Fe bonds show strong iron isotope dependence, which is a marker of an "intrinsic" isotope 
effect. This bulk iron isotope contribution to the local lattice dynamics suggests the local phonons play an 
important role in the (Ba, K)Fe 2 As 2 superconductor and a three-dimensional electron-phonon interaction 
should be considered to understand the superconducting mechanism of iron-based superconductors. 



Since the discovery of cuprate superconductors in 1986, one of the still open and challenging problems in 
condensed matter physics is the origin of the superconductivity in high-temperature superconductors 
(HTSCs) 1 . Two decades of efforts based on thousands of experimental and theoretical researches pointed 
out that electron-electron interaction (EEI) rather than the electron-phonon interaction (EPI) plays the key role 
to trigger the superconductivity in HTSCs materials 2 ' 3 . However, iron-based superconductors the new family of 
HTSCs discovered in 2008 does not show an agreement with this assumption 4 " 6 . Indeed, a clear iron isotope effect 
on T c was experimentally observed 7 " 10 . Chen et al. 7 firstly reported a large and positive iron isotope effect (Fe-IE) 
on the T c of both LnFeAs(0, F) and (Ba, K)Fe 2 As 2 superconductors. Khasanov et al. 10 found a larger Fe-IE 
exponent, well beyond the BCS value. Finally Shirage et al 8 ' 9 observed an inverse iron isotope effect in the (Ba, 
K)Fe 2 As 2 superconductor. Recently Yanagisawa et al. u presented a model explaining the sign inversion of the IE 
exponent. In his model, the (Ba, K)Fe 2 As 2 system is characterized by a multiband superconductivity originating 
by two different paring channels: a phononic contribution and an antiferromagnetic (AF) fluctuation mech- 
anism. However, his approach has been questioned by Bussmann-Holder and Keller 12 that claimed that an 
inversion of the Fe-IE exponent could not occur. More recently, Khasanov et al 13 tried solving the discrepancy 
in the Fe-IE measurements by decomposing the exponent in two terms: one associated to structural changes and 
one intrinsic. They also pointed out that the Fe isotope substitution may induce small structural changes which, in 
turn, may affect T c . Actually, the controversy still on points out that: 

1) lattice vibrations (i.e., phonons) play an important role in the superconductive mechanism(s) of iron-based 
superconductor materials; 

2) the understanding of the IE in these iron-based materials remains a challenging issue to clarify. 

Therefore, additional and more accurate experimental methods to measure the IE, probing not only mac- 
roscopic averaged properties such as those based on T c determinations, but also microscopic information of the 
system are required. 
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Here we present an investigation of the Fe-IE in the (Ba, K)Fe 2 As 2 
superconductor system by temperature- dependent extended x-ray 
absorption spectroscopy (T-EXAFS). Due to the correlation between 
the absorber and the backscatter through the emitted photoelectrons 
excited by x-rays, EXAFS probes the local structure around the 
selected photoabsorber, such as the number of neighboring atoms 
N, the average bond length between absorbing atoms and neighbor- 
ing atoms <r> and the variance a 2 1415 . The Debye- Waller factor a 2 
corresponds to the parallel mean square relative displacement 
(MSRD) among two atoms and it is the sum of two contributions: 

(i) 



i ^2 



The first term is temperature independent and associated to the static 
disorder, e.g., to the dispersion of bond lengths neglecting the vibra- 
tional dynamics. The second (dynamical) contribution is due to 
thermal vibrations and can be correlated to the vibrational spectrum 
of the system. The latter is usually well described by the Einstein 
model 16 

\T E ] 



# 1 u 
— — — coth 

2k{i T E 



2T 



(2) 



where fi is the reduced mass of the atom pair, T is the temperature, k 
is Boltzmann's constant and fi is the Plank's constant. In the Eq. 2, 
Te = (h/k)a)E is the Einstein temperature and co E is the correspond- 
ing Einstein frequency. The o vi ^ 2 contribution is different from zero 
also at 0 K. 

Measurements of T-EXAFS may allow, via the analysis of the 
Debye-Waller factor, to separate the static and the thermal disorder 
contributions present in a sample. The first term is relative to the 
local structure and affects the electronic and the magnetic behavior of 
the system, while the second depends only by the local vibration 
behavior (i.e., phonons). By using the T-EXAFS technique, we may 
evaluate at the microscopic scale the contribution of the Fe-IE on 
both the static structure and the lattice dynamics in iron-based SCs. 
Because EXAFS is an element- specific technique, in these systems we 
may independently achieve information of the individual atomic 
pairs such as the Fe-As and the Fe-Fe pairs. Powder samples can 
be also used in an EXAFS experiment and single crystals are not 
necessary to accurately determine the mean-square relative displace- 
ment (MSRD). 

Both Fe and As K-edge extended x-ray absorption fine structure 
spectra of the Bao.6Ko.4Fe 2 As 2 system by replacing 56 Fe ( 56 Fe-sample) 
with the isotope 54 Fe ( 54 Fe- sample) were measured vs. temperature 
down to 10 K. After an accurate analysis and comparison of EXAFS 
spectra of the two isotopic substitutions we point out that it does not 
induce detectable change of the average distances of both Fe-As and 
Fe-Fe atomic pairs while slightly modulates the Fe-As and Fe-Fe pair 
distributions. Moreover, it affects the local vibration of both Fe-As 
and Fe-Fe bonds. 

Results 

Fig. 1(a) compares As K-edge EXAFS data at 10 K and at room 
temperature (RT) showing a clear damping of the EXAFS oscillations 
vs. temperature. This behavior of the local structure is also clear 
looking at the comparison of the Fourier transforms (FTs) of the 
EXAFS spectra in Fig. 2(b). In this iron-based superconductor with 
the 122 -structure no anomalies are detected at low temperatures, 
especially around T c . In addition to our investigation, to our know- 
ledge other two manuscripts published temperature-dependent 
EXAFS measurements on iron-based superconductors. They are 
from the Oyanagi and Tyson teams. Both investigated the (1111) 
iron-based superconductors LaFeAsO!_ x F x obtaining different 
results. Oyanagi et al. found an anomalous behavior across T SDW 
and T c 17 . No such fine behavior has been observed in the data of 
Tyson 18 . We are convinced that in addition to always possible 
systematic errors or limited temperature points, the differences 



observed can be associated to the different F concentration: x = 
0.07 and x = 0.11. Actually, the temperature- dependence of Fe-As 
and Fe-Fe MSRD of our system is similar to that observed by Tyson et 
al, but our work is investigating a different type of iron-based super- 
conductors: the (Ba, K)Fe 2 As 2 system belonging to the (122) family 
whose structure is well different from pnictides of the (1111) family. 
On the contrary, as we will discuss in the next, the iron isotope effect 
can be clearly seen looking at the magnitude of the FTs, which is 
correlated with the a 2 term, i.e., the EXAFS Debye-Waller (DW) 
factor. 

In Fig. 1(c) and (d) we compare the first-shell peak (As-Fe pair) of 
the FTs of the 54 Fe- and 56 Fe-samples at 10 K and at RT, respectively. 
The position of the first-shell peak, i.e., the As-Fe pair in the FTs, is 
constant in the two samples and the same occurs for the Fe-Fe pair 
(panel d). Data point out that the Fe-IE does not affect the average 
distances of both Fe-As and Fe-Fe atomic pairs in the time-domain of 
the EXAFS experiments. A difference in the FTs associated to the 
isotopic substitution appears only looking at the magnitude of the 
first- shell peak. As it is well know, for a given crystal structure the 
magnitude of the FT peak depends by the a 2 term. As shown in 
Fig. 1 (c) a difference in the FT magnitude for the As-Fe pair is evident 
at the lowest temperature (10 K - panel c). It gradually reduces with 
the temperature and at RT the two curves almost overlap (panel d). 
To obtain a quantitative dependence of a 2 vs. T, we extracted the DW 
of both As-Fe and Fe-Fe pairs by fitting EXAFS signals vs. T from 10 
K to RT at both As and Fe K-edges, respectively. For the As K-edge, 
we considered the EXAFS signals of the nearest As-Fe bond. For the 
Fe K-edge, both nearest neighboring (Fe-As) and next-nearest neigh- 
boring (Fe-Fe) were filtered out. We have to underline here that the 
As-Fe correlation is identical to the Fe-As correlation in an EXAFS 
process. Therefore, in order to obtain exact information on the Fe-Fe 
bond we fix the distance and the DW factor of the Fe-As bond using 
parameters obtained by the As K-edge EXAFS fit during the analysis 
at the Fe K-edge EXAFS. Fig. 2a shows the quality of the EXAFS fit of 
the As K-edge EXAFS oscillations at RT. In Fig. 2b we compare the 
average distances of both Fe-As and Fe-Fe bonds of 54 Fe- and 56 Fe- 
samples, respectively, in the whole temperature range (10 ~ 300 K). 
Data are consistent with what showed in Fig lc and Id. We may 
emphasize that the result is different from the report of Khasanov 
et at. 10 on the FeSe!_ x system that investigated the temperature 
dependence of the lattice structure by Neutron Powder Diffraction 
(NPD). This work also shows a large Fe-IE (oc Fe = 0.81) by means of 
SQUID measurements. As claimed by Khasanov et al, it is the tiny 
modulation of the structural parameters caused by the isotope sub- 
stitution that induces the large Fe-IE, beyond the framework of the 
BCS theory. However, these findings indicate that a Fe isotope effect 
on T c is really present and related with the atomic or lattice structure 
in the Fe-based superconductors of the (11) family. In the super- 
conductor (Ba, K)Fe 2 As 2 [(122) family] we detected an Fe-IE only on 
the atomic distributions and not on bond distances. Moreover, data 
points out that the isotope substitution has no effect on lattice dis- 
tances in these samples of the (122) family. As a consequence, it 
represents an ideal model to characterize the isotope effect on the 
vibrational contributions, a fundamental issue to understand the 
electron-phonon interaction in a superconductor. 

In Fig. 3 we also show the dependence of the DW factors of the As- 
Fe and the Fe-Fe bonds vs. temperature for both the 54 Fe- and 56 Fe- 
samples. A clear "isotopic" contribution appears at low temperature 
while at high temperature, the dependence of the DW factor 
approaches the classical linear behavior. 

As outlined above, the a 2 term contains two contributions: the 
static disorder and the thermal disorder. In order to evaluate the 
contribution of the static disorder and of the vibrational frequency 
of the As-Fe and Fe-Fe bonds of isotopic substituted samples, 
we used the correlation Einstein model to fit the a 2 function vs. 
temperature: 
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Figure 1 | Temperature-dependent EXAFS data for the Fe isotope substituted samples, (a) Representative experimental EXAFS oscillations at the As 
(upper panel) and Fe (lower panel) K-edges of the 54 Fe-sample. (b) FTs of the As K-edge EXAFS data measured cooling the sample down to 10 K. 
(c) Comparison between the first-shell peak of the FTs of the 54 Fe- and 56 Fe-samples at 10 K. (d) The same of (c) at RT. 
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In Eq. 3 o\ is the static contribution to the bond and the quality of our 
simulations is shown by solid lines in the same figure. The corres- 
ponding values of the static contribution and of the vibrational fre- 
quency of the two samples are summarized in table 1 for both 
As-Fe and Fe-Fe bonds. Data show that the isotopic substitution 
gives rise to a weak modulation of the distribution of the interatomic 
distances of both Fe-As and Fe-Fe pairs. Moreover, the Fe-IE affects 
slightly more the static structural disorder of the Fe-Fe pair (it 
changes from 0.0056 to 0.0053 A 2 ) than that of the As-Fe (or Fe- 
As) pair that decreases from 0.0025 to 0.0023 A 2 . In addition, the 
static disorder of the Fe-Fe pair is much larger than that of the Fe-As 
pair. Bonds with a larger o\ value (i.e., the Fe-Fe pair) and modu- 
lation due to the isotopic substitution, point out that these bonds are 
more flexible. This result is in agreement with other experimental 
data 19 and because decreasing temperature down to the structural 
transition temperature (Ts) the Fe-Fe bond will undergo to a split- 
ting 20 , data are also consistent with the occurrence of a structural 
phase transition in the iron layer of the BaFe 2 As 2 system. Also by 
using T-EXAFS technique, in the K 0 .8Fei. 6 Se 2 [122-chalcogenide] a 
similar lattice property has been detected. Also this work points out 
that the Fe-Fe bond is more flexible than the Fe-As bond 21 . The iron 



isotope effect we observe on o\ is due to the static structural modu- 
lation and corresponds to electronic and magnetic properties of the 
iron-based SCs. 

A second important result we obtained by the analysis of the T- 
EXAFS experiments is that the Einstein temperature (T E ) that is a 
probe of the vibrational frequency of the system, exhibits a shift for 
both the As-Fe and Fe-Fe bonds induced by the iron isotope substi- 
tution. The T E of the As-Fe bond changes about 1.3% from 331.5 to 
327.1 K after the substitution with 54 Fe in the Bao.6Ko.4Fe 2 As 2 . In the 
Fe-Fe bond, T E shifts -1.5% from 169.0 to 166.5 K. Both Fe-Fe and 
Fe-As vibrations have then a similar, but not negligible contribution 
to the superconductive mechanism of iron -based SCs. Although the 
detected changes are small, similar changes of T E occur for both Fe- 
Fe and Fe-As bonds. Actually, the Einstein temperature T E can be 
associated to an effective force constant k 0 normalized to the reduced 
mass |i: 

T E = ^kJji. (4) 

Actually, because the Einstein temperature is lower for the Fe-Fe 
bond, the latter is supposed to be softer than the Fe-As one. This 
experimental result points out a greater flexibility of the iron layers in 
the (Ba, K)Fe 2 As 2 system. To evaluate the contribution of the isotope 
effect on the local vibration of the Fe-Fe and As-Fe correlations, we 
introduced the coefficient a 
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Figure 2 | EXAFS fits and result of the bond distances vs. T. 

(a) Experimental (black circles) and fit (red lines) EXAFS oscillations of 
As-Fe (upper curve) and Fe-Fe (lower panel) shells, respectively, for the 
54Fe-sample at RT. (b) Temperature dependence of Fe-As bond distances 
and Fe-Fe bond distances for the 54 Fe- and the u Fe-samples, respectively. 
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Figure 3 | DW factors a 2 for the bond of the 54 Fe- and the 56 Fe-samples 
vs. T. (a) As-Fe bond; (b) Fe-Fe bond. Symbols are experimental data 
extracted from the EXAFS fits while solid lines are simulations based on the 
Einstein model. Bars refer to statistical errors. 



very similar to the definition of the T c isotope effect coefficient 22 . 
Here T E is the Einstein temperature obtained by the EXAFS experi- 
ments, regarding both Fe-Fe and Fe-As bonds (Table 1) while M is 
the atomic mass of the isotopic substitution. Using the coefficient a 
we may obtain an independent evaluation of the iron-isotope- effect 
for both As-Fe and Fe-Fe bonds. Using values listed in table 1 we got 
values of 0.37 and 0.41 for the As-Fe and Fe-Fe bonds, respectively, 
very close to the T c iron-isotope-effect coefficient of 0.37 reported 
by the Chen's group 7 . Our analysis supports the hypothesis that 
the electron-phonon interaction plays a role in the supercon- 
ductive mechanism of iron-based superconductors and that both 
Fe-Fe and Fe-As vibrations contribute to the superconducting 
mechanism. 

Because Fe-Fe bonds lie in the superconducting layers (a-b plane) 
while Fe-As bonds are oriented along the c axis, our result confirms 
in the iron-based superconductors the presence of a weak bulk 
iron isotope effect. In addition, it points out the presence in iron- 
based superconductors of a not negligible 3D electron-phonon inter- 
action as shown in Fig. 4. This result can be also used to interpret 
the nearly isotropic superconducting properties of iron-based 
superconductors 23 , systems well different from layered cuprate 
superconductors. 



Discussion 

An accurate EXAFS analysis of the (Ba, K)Fe 2 As 2 system vs. temper- 
ature shows that an iron isotope effect correlated to the static struc- 
tural distribution and with local vibrations can be detected. Changes 
of the static structural distribution may occur via changes of the 
inter- atomic distances and of the atomic distribution (or of the dis- 
order). The Fe-IE does not induce appreciable changes in the inter- 
atomic distances while the atomic distribution is slightly affected. At 
present, it remains however extremely difficult identify a clear cor- 
relation among local disorder contributions (present in a powder 
sample) and superconducting mechanisms. In addition, we clearly 
showed the existence of a Fe-IE contribution to local vibrational term 
for both Fe-As and Fe-Fe bonds. Their isotopic exponents are in 
agreement with previous experiments 7 looking at the iron isotope 



Table 1 | Values of local static distortions and Einstein tempera- 
tures associated to the As-Fe and Fe-Fe bonds in the 54 Fe- and 
56 Fe-samples 

bonds Isotope o 0 2 {A 2 ) T E (K) a 

AsTe ^Fe 0.0025 ± 0.0002 331.5 ± 2.0 0.37 

56 Fe 0.0024 ± 0.0002 327.1 ± 2.0 

Fe-Fe 54 Fe 0.0056 ± 0.0002 1 69.0 ± 1 .0 0.41 

56 Fe 0.0053 ± 0.0002 1 66.5 ± 1 .0 
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Figure 4 | Layout of Fe-Fe and Fe-As vibrations in the iron-based superconductors. Fe-Fe bonds are in the superconducting layers and their vibrations 
attract electrons in the c direction, while Fe-As bonds are oriented along the c axis and their vibrations affect electrons in the a-b layers. 



effect in these systems. This result represents the first experimental 
link between the Fe-IE of Tc experiments and the superconducting 
mechanism at a microscopic level. 



Methods 

Isotopically substituted polycrystalline samples were synthesized by the solid state 
reaction method described in detail in Ref. 7. Previous investigation of this super- 
conductor showed a large iron isotope effect on the T c 7 . As and Fe K-edge absorption 
spectra were collected using a Si (111) double crystal monochromator at the XAFS 
station of the 1W1B beamline of the Beijing Synchrotron Radiation Facility (BSRF). 
The storage ring was working at 2.5 GeV with an electron current decreasing from 240 
to 160 mA within about 8 hours. Samples were ground into fine powders and brushed 
onto tapes that were stacked together to give approximately one x-ray-absorption 
lengths at their corresponding metal edges. Samples were cooled by a closed-cycle 
helium refrigerator with an error bar on the temperature determination of about 0.1 
K. Data were processed and analyzed using both UWXAFS 24 and FEFF 8.0 
packages 25,26 to generate scattering amplitudes and phase shifts starting from available 
XRD refinements (see supplementary materials of Ref. 7). 
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